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t-h ' Abstract 

> : 

QO , We show that the problem of designing a quantum information error correcting pro- 

cedure can be cast as a bi-convex optimization problem, iterating between encoding and 
recovery, each being a semidefinite program. For a given encoding operator the problem 

O i s convex in the recovery operator. For a given method of recovery, the problem is con- 

vex in the encoding scheme. This allows us to derive new codes that are locally optimal. 
We present examples of such codes that can handle errors which are too strong for codes 

^ I derived by analogy to classical error correction techniques. 

1 Introduction 

P . 
cr 

Quantum error correction is essential for the scale-up of quantum information devices. A the- 
ory of quantum error correcting codes has been developed, in analogy to classical coding for 
noisy channels, e.g., IIS95 1 IS961 lG%l IKL97I . Recently IRW051 and IYHT05I did this by posing 
error correction design as an optimization problem with the design variables being the process 
matrices associated with the encoding and/or recovery channels. Using fidelity measures leads 
naturally to a convex optimization problem, specifically a semidefinite program (SDP) IB V04H . 
The advantage of this approach is that noisy channels which do not satisfy the standard assump- 
tions for perfect correction can be optimized for the best possible encoding and/or recovery. 

In [RW05| the power-iteration method was used to find optimal codes for various noisy 
channels, by alternately optimizing the encoding and recovery channels. In contrast, here we 
apply convex optimization via SDR and similarly iterate between encoding and recovery. For a 
given encoding operator the problem is convex in the recovery. For a given method of recovery, 
the problem is convex in the encoding. We further make use of Lagrange Duality to alleviate 
some of the computational burden associated with solving the SDP for the process matrices. 
The SDP formalism also allows for a robust design by enumerating constraints associated with 
different error models. We illustrate the approach with an example where the error system does 
not assume independent channels. 
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An intriguing prospect is to integrate the results found here within a complete "black-box" 
error correction scheme, that takes quantum state (or process) tomography as input and iterates 
until it finds an optimal error correcting encoding and recovery. 



2 Quantum Error Correction 
2.1 Standard model 

A standard model [NC00 ( §10.3] of an error correction system as shown in the block diagram 
of Figure[T]is composed of three quantum operations: encoding C, error £ , and recovery 1Z. 
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Figure 1 : Standard encoding-error-recovery model of an error correction system. 

The input, p s , is the n s x n s dimensional density matrix which contains the quantum infor- 
mation of interest and which is to be processed. We will refer to ps as the system state or the 
unencoded state. The output of the encoding operation is pc, the nc x nc dimensional encoded 
state. The error operator, which is also the source of decoherence, corrupts the encoded state 
and returns ac, the nc x nc "noisy" encoded state. Finally, pn is the tir x tir dimensional 
recovered state. The objective considered here is to design (C, 1Z) so that the map ps — > Pr is 
as close as possible to a desired n s x n s unitary L s . Hence, pr has the same dimension as ps, 
that is, n R = n s . For emphasis we will replace p R with p s . 

Although it is possible for £ to be non-trace preserving, in the model considered here, all 
three quantum operations are each characterized by a trace-preserving operator-sum-representation 
(OSR): 

pc = C(p s ) = CcPsCl C & = J ^> C * e Cn ° XnS 



a C = £( Pc ) = J2 E eP0El J2 E e E e = I no S e £ 



(1) 



%sxnc 



PS = K{a C ) = Yl R r°cRl R r R r = 7 »C> R r G 

r r 

These engender a single trace-preserving quantum operation, S, mapping p s to ps, 

Ps = S(p s ) = Y S rec p s S}. ec 

r,e,c 

(2) 

Srec = RrE e C c 6 C S S =?■ ^ ^ S]. ec S rec = I ris 

r,e,c 

Before we describe our design approach we make a few remarks about the error source and 
implementation of the encoding and recovery operations. 



2 



2.2 Implementation 

Any OSR can be equivalently expressed, and consequently physically implemented, as a uni- 
tary with ancilla states [NCOO, §8.23]. An equivalent system-ancilla-bath representation of the 
standard error correction model of Figure[T]is shown in the block diagram of FigureEl 
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Figure 2: System-ancilla-bath representation of standard encoding-error-recovery model of er- 
ror correction system. 

For the encoding operator, C, the encoding ancilla state, \0ca), has dimension uca, and 
hence, the resulting encoded space has dimension n c = n s n C A- The encoding operation is 
determined by U c , the n c x n c unitary encoding operator which produces the encoded state 
<>c = U c (ps®\0 C a)(Oca\)uI . 

For the error system, E, the ancilla states are engendered by interaction with the environ- 
ment, or the term used here, the bath. Here we will ignore complications associated with an 
infinite dimensional bath. The error system is thus equivalent to the he X ue unitary error op- 
erator Ue with uncoupled inputs, pc the encoded state, and p R , the % x tib bath state. Thus, 
n>E = nsncAns- The noisy encoded state ac, is the n c x n c reduced state obtained by tracing 
out the bath from the output of U E , that is, o c = Tr B U E {pc ® Pb)U e . 

The recovery system 1Z has additional ancilla \0ra) of dimension Ura- Uh is the tir x tir 
unitary recovery operator with n R = nsncA^RA and with a R the n R x n R full output state 
vr = Ur(ctc ® \®ra) (Ora\)U r . The n s x n s reduced output state, ps, is given by the partial 
trace over all the ancillas, the bath having been traced out in the previous step. Specifically, 
ps = Tr A a R 



Caveat emptor The "real" error correction system is unlikely to be accurately represented 
by the system shown in Figure El but rather by a full system-ancilla-bath interaction [ALZ05|. 
As shown in the block diagram in Figure|3J Uqab is the uqab x uqab unitary system-ancilla- 
bath operator, IOcaO^a) is the total ancilla state of dimension uca^ra and ps is the bath 
state. The reduced system output state, ps, is obtained from the full output state pqab by 
tracing simultaneously over all the ancilla and the bath, p s = Tr ab Pqab- At this level of 
representation, there is no distinction between the n C A ancilla states used for encoding and the 
n RA ancilla states used for recovery. The internal design, however, may be constructed with 
such a distinction. 
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Figure 3: System-ancilla-bath representation of error correction system. 

2.3 Optimal error correction: maximizing fidelity 

As stated, the goal is to make the operation p s — > p s be as close as possible to a desired unitary 
operation L s . Measures to compare two quantum channels are typically based on fidelity or 
distance, e.g., HGLN05I . HKGBR06I . Let S denote a trace-preserving quantum channel mapping 
n-diniensional states to n-dimensional states, 



(3) 



The following fidelity inequalities hold: 

/mixed ^ /pure ^ /avg 

where 

/mixed = mm p ^ k |Tr S k p\ 



(4) 



pure 



min l^) E fc KV>l<W)l 



(5) 



/avg — Efc |Tr5fc| 



All are in [0,1] and equal to one if and only if S(p) = p. From ©, S = L^TZSC with OSR 
elements S k = S rec = L%R r E e C c e C nsXns . Thus / avg = 1 <^> S = l s <^> S = TZ£C = L s . 

Given S, not all these fidelity measures are easy to calculate. Specifically, / m i xe d is a convex 
optimization over all densities, that is, over all p 6 c nxn , p > 0, Tr p = 1, and hence, can 
be numerically obtained. Calculation of / avg is direct. Calculating / pure is, unfortunately, not a 
convex optimization over all pure states If, however, the density associated with / m i xe d is 
nearly rank one, then / mixcd w / pure . 

As a practical matter, when dealing with small channel errors, it does not matter which 
fidelity measure is used. Therefore it is convenient to use / avg , as it is already in a form explicitly 
dependent only on the OSR elements. In [RW05 1 / avg was also used as the design measure, but 
specific convex optimization algorithms were not proposed. In [ YHT05 1 a similar optimization 
was proposed using a distance measure to obtain the recovery given the encoding. 

We now focus on the optimization problem, 



maximize f avg (TZ,C) 
subject to J2 r R\R r 



(6) 



L nc i 



E, c\c c = i, 



I'S 
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The optimization variables are the OSR elements {C c } and {R r }- As posed this is a difficult 
optimization problem. The objective function is not a convex function of either of the design 
variables. In addition, the equality constraints are quadratic, and hence, not convex sets. The 
problem, however, can be approximated using convex relaxation, where each nonconvex con- 
straint is replaced with a less restrictive convex constraint [BVOfl. This finally results in a 
bi-convex optimization problem in the encoding and recovery operator elements which can be 
iterated to yield a local optimum. As we show next, iterating between the two problems is guar- 
anteed to increase fidelity of each of the relaxed problems. When the iterations converge, all 
that can be said is that a local solution has been found. 

3 Optimal error correction via bi-convex relaxation 

3.1 Process matrix problem formulation 

Following the procedure used in quantum process tomography HNC001 §8.4.2], [KWR04J we 
expand each of the OSR elements R r G Qn s xn c anc j e Qn c xn s m a set f basis matrices, 
respectively, for C nsXnc and C ncXns , that is, 

= J2i X ci Bci 

where { B m G C nsXnc , B Ci G C n ° xns \ % = 1, . . . , n s n c } and the {x ri } and {x ci } are com- 
plex scalars. Problem © can then be equivalently expressed as, 

maximize f avg (TZ, C) = f avg (X R , X c ) = J2 ijki (Xr)^ (X c )h F ijM 
subject to ( x R)ij B Ri B Rj = In c 
J2u ( x c)k£ Bc k B Ci = I ns 

(8) 

( X iij ij %ri%rj 

{X c )u = E c XckXrf 

Fijke = E e (Tr L\B m E e B Ck ){Tr L\B Rj E e B Ci y /n| 

The optimization variables are the process matrices X R , X c G C nsncXnsnc and the scalars 
{x r i} and {x C i}. The problem data which describes the desired unitary and error system is con- 
tained in the {F ijke }. The equality constraints {Xr)^ = E r x nKj and {Xc)u = E c ^cfc^ 
are both quadratic, exposing again that this is not a convex optimization problem. We do not 
explore the possible simplifications that can occur in these expressions if the basis matrices are 
chosen prudently, e.g., Tr BjBj = 5^. 

3.2 Design of 1Z given C and £ 

In this section and in the remainder of the paper we set the desired logical operation to identity, 
i.e., L s = I s ', just error correction not correction and computation. This is without loss of 
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generality as a desired logical operation can be added everywhere. 

Suppose the encoding C is given (and L s = I s ). Then optimizing only over 1Z in © can be 
equivalently expressed as, 

maximize f avg {n, C) = f avg {X R , C) = Tr X R W R (S, C) 
subject to ( x R)ij B m B Rj = l n c 

(W R (£,C)h = E ck e F m = £ ec (Tr B Ri E e C c )(Tr B Rj E e C c )* jn\ 

(9) 

The optimization variables are the matrix X R E C nsncXnsnc and the scalars {x ri }. The prob- 
lem data is contained in the positive semidefinite matrix W R (£,C) E Q n sncxn S n c _ 
jective function is now linear in X R , which is of course a convex function. However, each of 
the equality constraints, (X R )ij = J2 r x ri x*j is quadratic, and thus does not form a convex 
set. This set of quadratic equality constraints can be relaxed to the matrix inequality constraint, 
X R > 0, that is, X R is positive semidefinite, a convex set in the elements of X R . A convex 
relaxation of © is then, 

maximize Tr X R W R (£, C) 

(10) 

subject to X R > 0, ( X R)ij B m B Rj = 4 C 

This class of convex optimization problems is referred to as an SDP, for semidefinite program 
HBV04I . 1 For a given encoding C, the optimal solution to the relaxed problem (TlOb . X^ x " opt , 
provides an upper bound on the average fidelity objective in © or ©. From the fidelity in- 
equalities ©, we can derive a lower bound. Specifically, the (unknown, possibly unknowable) 
solution to the original problem ©, is bounded as follows: 

/mixed (7^ lx - op \ C) < max / pure (ft, C) < /av g (^ rlx - op \ C) = Tr X^WrS^, C) (1 1) 

where ^ rlx - opt is the OSR with elements {i^ lx - opt } obtained from Xjj x " opt via the singular value 
decomposition, 

nsnc 

x Hx_opt = y S y f ^ R ri x _ op t = ^ V ir B Ri , r = 1, . . . , n s n c (12) 

i=i 

where V E C nsncXnsnc is unitary and S = diag(si ■ • • s nsnc ) with singular values in decreas- 
ing order, s 1 > s 2 > ■■ ■ > > 0. 

'A standard SDP is to minimize a iinear objective function subject to convex inequalities and linear equalities. 
The objective function in (1101 is the maximization of a linear function which is equivalent to the minimization of 
its' negative, and hence, is a linear objective function. 
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3.3 Design of C given 7Z and 8 

Repeating the previous steps, optimizing only over C in © can be equivalently expressed as, 
maximize / avg (ft, C) = / avg (X c , K) = Tr X C W C (S, K) 

subject to J2k,i ( x c)ke B Ck B ce = I ns 

£ e r (Tr B Ck R r E e )(Tr B C eR r E e )*/n s 

(13) 

The optimization variables are the matrix X c E C™ sncXrisnc and the scalars {x ci } with all the 
problem data contained in the symmetric positive semidefinite matrix Wc (S,TZ) E C nsnc xn s n c 
In this case, however, the basis matrices, are n c x ns- Repeating the previous procedure 

of relaxing the quadratic equality constrain to Xc > 0, we obtain the convex relaxation of ® 
as the SDP, 

maximize Tr X C W C (S,TZ) 

(14) 

subject to X c >0, Eij ( x ch B^Bcj = I ns 

Analogously to (fT2t . for a given recovery 1Z, the (unknown, possibly unknowable) solution to 
the original problem ©, is bounded as follows: 

/mixed (n, c rlx - opt ) < max f pum (n, c) < f avg (n, c rlx - opt ) = tv x c x - opt w c (£, n) (15) 

where C rlx - opt is the OSR with elements {CJ lx - opt } obtained from X c x - opt via the singular value 
decomposition, 

nsnc 

x rix.o P t = vsy] ^ ^ix.opt = ^ V ic B Ci , c = 1, . . . , n s n c (16) 

i=i 

where V E C nsncXnsnc is unitary and S = diag(sx ■ • • s nsnc ) with singular values in decreas- 
ing order, s 1 > s 2 > ■■ ■ > s nsnc > 0. 

3.4 Iterative bi-convex algorithm 

Proceeding analogously as in [RW05 1, the two separate optimizations for C and 1Z can be com- 
bined into the following iteration. 

initialize encoding C and stopping level e 
repeat 

1 . optimize recovery 
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(a) compute X R as solution to: 

maximize Tr X R W R (£, C) 

subject to X R > 0, ( x n)ij B m B Rj = l n c 

(b) use (fT2ll to compute 1Z* from X R 

2. optimize encoding 

(a) compute X* c as solution to: 

maximize Tr X C W C (S, K*) 

subject to X R >0, J2i,j( x ch B a B cj = 4 S 

(b) use (fToT) to compute C* from X£ 

3. compute change in fidelity 

4. reset 

C = C* 

until 

The algorithm returns (7?*, C*). The optimization in each of the steps is a convex optimization 
and hence fidelity will increase in each step, thereby converging to a local solution of the joint 
relaxed problem. This solution is not necessarily a local solution to the original problem © or 
®. However, the upper and lower bounds obtained will apply. The optimization steps can be 
reversed by starting with an initial recovery and then starting the iteration by optimizing over 
encoding. 



3.5 Decoherence resistant encoding 



If the sole purpose of encoding is to sustain the information state ps, then the desired operation 
is the identity (Ls = Is) and the recovery operation in Figure [T]is simply the partial trace over 
the encoding ancilla, that is, 



ps = H{crc) = Tt C a oc 



Tr(ac)[i,i] 
Tr (o"c)[ ns ,i] 



Tr Oc)[i,n s ] 



Tr (<T )[n 



S,ns\ 



(17) 



where the (ac)u,j] are the n 2 s sub-block matrices of <Tc, each being uca x "^ca- Hence, the OSR 
elements of 1Z are given by 



(Rr) 



1 j = (i- l)n C A + r 
else 



ncA, j = l,...,ns 



(18) 
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For a given error £, finding an optimal encoding by solving (TRb is equivalent to finding a 
decoherence-resistant-subspace. If there is perfect recovery, then we have found a decoherence- 
free-subspace HLCW9 8I. In [ZL04], this problem was considered using / purc , the pure state 
fidelity. 

3.6 Robust error correction 

The bi-convex optimization can be extended to the case where the error system is one of a 
number of possible error systems, that is, 

£e{£ a \a = l,...,£] (19) 

where each £ a has OSR elements {E ae }. The worst-case fidelity design problem, by analogy 
with ©, is then: 

maximize mine, f avs (7Z,£ a ,C) = Y]^^ (X R )ij (Xc) kl Faijkl 

subject to X R , X c constrained as in dU) (20) 

Faijke = J2 e LsB Ri E ae B C k)(Tr LgB Rj E ae B C e)* /n 2 s 

Iterating as before between 1Z and C results again in separate convex optimization problems, 
each of which is an SDR Specifically, for a given encoding C, a robust recovery is obtained 
from, 

maximize min a TV X R W R (£ a , C) 

(21) 

subject to X c > 0, ( x c)i 3 B Ci B Cj = I ns 

Similarly, for a given recovery 1Z, a robust encoding is obtained from, 

maximize min a Tr XcWc(£ a , Tty 
subject to X R > 0, £\ j {X R )ij B Ri B Rj = I nc 

4 Computing the solution: Lagrange Duality 

The main difficulty with embedding the OSR elements into either X R or X c is scaling with 
qubits. Specifically, the number of design parameters needed to determine either Xc or X R 
scales exponentially with the number of qubits. Although exponential scaling at the moment 
seems unavoidable, we show in this section that solving the dual SDPs associated with either 
dTOb or dT~4T> requires many fewer parameters, and thus engenders a reduced computational bur- 
den. 

The convex optimization problems dTOb and (fl4ll are both SDPs of the form, 

maximize Tr XW 

subject to X > 0, XijB\ B i = J ™ 



(22) 
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with optimization variable X = E C nxn , n = rm for some integer r, and with each basis 
matrix Bi E C" xm , We will refer to this as the primal problem. From IBV041 §11.8.3], for 
the standard SDP: minimize c T x, subject to F + XiFi > with x E R p and Fi = Ff E 
R 9 * 9 , the computational complexity using a primal-dual algorithm is max {pq 3 ,p 2 q 2 ,p 3 } flops 
(floating point operations) per iteration step where typically 10-100 steps are required in the 
algorithm. Accounting for the linear (matrix) equality constraint and the Hermiticity of X, the 
number of real optimization variables in (123)1 is p = n 2 — m 2 = (r 2 — l)m. The dimension 
of the linear matrix inequality is q = n = rm. This gives the computational complexity as 
p 2 q 2 = r 2 (r 2 — l) 2 m 6 flops per iteration. 

Solving (ITOb for X R , gives n = n s n c , r = n s , m = n c = n s n C A for n s( n s ~ l) 2n cM 
flops. Solving (IT4T> for X c , gives n = n s n c , r = n c , m = n s for n s s {n 2 s n 2 CA — l) 2 n 2 CA flops. 
Exponential growth in computation occurs becasue each of these dimensions are exponential in 
the number of qubits, i.e., , ns = 2 qs , nc = 2 qs+qcA , and so on. 

The computational burden can be somewhat alleviated by appealing to Lagrange Duality 
Theory [BV04, Ch.5] which provides a means for establishing a lower bound on the optimal 
objective value, establishing conditions of optimality, and providing, in some cases, and this 
case in particular, a more efficient means to numerically solve the original problem. In Ap- 
pendixlAlwe show that the dual problem associated with the primal problem (l23l is, 

minimize Tr Y 



subject to K(Y) - W > 0, K^Y) = Tr B\B^Y 



(24) 



with optimization variable Y = Y^ E C mxm . The number of (real) optimization variables for 
the dual problem is then at most m 2 . The dual problem is also an SDP and from the previous 
formula therefore requires r 2 m 6 flops per iteration, a reduction in flops per iteration from the 
primal by a factor of (r 2 — l) 2 . We show in Appendix lAl that if (X opt , Y opt ) solve the primal 
and dual problems respectively, then: 

TrX opt W = TrF opt 
(K(Y opt ) - W)X opt = ( } 

The second equation above together with the linear equality constraint in (|2"31 can be used to 
obtain the primal solution X opt from the dual solution Y opt . That is, solve for X opt from the set 
of linear equations, 

(K(Y opt ) — W)X opt = 

E^X^BjB, = I m (26) 

Solving this type of linear set of equations is an eigenvalue problem and thus requires on the 
order of no more than n 2 flops IIGL83I. Thus the dual takes r 2 m 6 flops per iteration plus 
r 2 m 2 flops one time to convert from dual to primal. This is in comparison to the much larger 
r 2 (r 2 — l) 2 m 6 flops per iteration for the primal alone. Neglecting the dual to primal conversion, 
the speed-up in flops per iteration to calculate Xr is approximately [n 2 s — l) 2 , n s = 2 qs and 
for X c it is (n 2 c - l) 2 , n c = 2 qs+qcA . 



5 Example 

In this illustrative example, the goal is to preserve a single information qubit using a single 
ancilla qubit. Thus, the desired logical gate is the identity, that is, L s = I 2 , with n s = ncA = 2, 
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and hence, n c = 4. We made two error systems, £ a and £ b , by randomly selecting the unitary 
bath representation as shown in Figure |3 as follows: Each error system has a single qubit bath 



m E xn E 



n E 



state, |0) B , thus % = 2. The Hamiltonian for each system, He = H E e C 
n c n B = 8, was chosen randomly and then adjusted to have the magnitude (maximum singular 
value) H-HbH — 5e — 0.75. Then the unitary representing the error system was set to Ue = 
exp(— iHe) and from this the corresponding OSR £ was computed. The OSR elements to three 
decimal points for the two random systems are as follows. 



£ a < 



E a l 



E 



a2 



Ebi 



Eb2 



0.9 - 0.049?' 

-0.159 + 0.148? 
0.167 + 0.061? 
0.124 + 0.137* 

0.053 - 0.063? 
-0.168 - O.Oli 
0.119 + 0.053?' 
0.123 + 0.066? 

0.943 + 0.018? 

0.107 + 0.062? 
-0.025 - 0.042? 
-0.017-0.094? 

0.07?' 
-0.003 -0.147?' 
0.049 + 0.084? 
-0.191 + 0.095?' 



0.193 + 0.194? 
0.887 - 0.046? 
-0.07 + 0.004? 
0.167-0.155? 

-0.034 + 0.082?' 
0.141 +0.073? 
0.207 - 0.02? 
0.027 + 0.008? 

-0.14-0.024? 
0.876 + 0.068? 
0.122 + 0.073?' 
0.095 + 0.035? 

-0.2 - 0.082? 

0.138-0.155? 
-0.149 + 0.217?' 
-0.081 - 0.097?' 



-0.161 + 0.039?' 

0.148 -0.025? 

0.905 + 0.161? 
-0.203 + 0.118? 

0.148 -0.085? 
0.008 + 0.091? 
0.043 + 0.01? 
0.07 - 0.058?' 



0.076 - 
-0.06 - 
0.889 - 
-0.032 



0.081?' 
■ 0.021? 

0.035? 
- 0.089? 



0.028 - 0.083?' 
0.202 + 0.306? 
0.143 -0.04?' 
0.007 - 0.127? 



-0.135 + 0.156?' 
-0.168- 0.081? 

0.16 + 0.125? 

0.844 - 0.26? 

0.13 - 0.076?' 
-0.074 - 0.024? 
-0.063 -0.21? 

0.098 - 0.11?' 

0.04 - 0.163? 
-0.127 + 0.06?' 
0.043 - 0.078?' 

0.88 + 0.113? 

0.179 + 0.206?' 
0.045 - 0.134? 
0.024 + 0.174? 
0.035 - 0.167? 



Neither of these error systems is of the standard type, e.g., there is no independent channel 
structure. The choice of 5e = 0.75 is perhaps extreme, but is motivated here by our desire 
to demonstrate that the optimization procedure can handle errors that are beyond the range of 
classically-inspired quantum error correction. For this particular set of error systems, we do not 
know if there exists an encoding/recovery pair limited to using a single encoding ancilla state 
which can bring perfect correction. This also motivates the search for the still elusive black-box 
error correction discussed in the introduction. 

For each of the error systems we ran the bi-convex iteration 100 times starting with the 
initial recovery operator 1Z given by the partial trace operation (fTSl . Denote (H a i,C a i) and 
(T^aioO) C a ioo) as the 1st and 100th iteration pairs optimized for £ a , and similarly (7£&i, C b i) and 
(T^-feiocb C&ioo) as the 1st and 100th iteration pairs optimized for £ b . Table[T]shows the average 
fidelities fa, vg (1Z, £, C) for some of the possible combinations. 

As Table [T] clearly shows, fidelity tuned for a specific error, either £ a or £ b in this example, 
saturated to the levels shown (0.9997) in about 100 iterations. However, neither of the optimized 
codes are robust. Each does very poorly when the error is different then what was expected. By 
raising the number of ancilla it is of course possible to make the system robust. This, however, 
introduces considerable complexity. What the table suggests is that an alternate route is to tune 
for maximal fidelity, say, in a particular module. This of course can only be done on the actual 
system. 



11 



Type 


n, C 


£a 


£ b 


Optimal encoding, no recovery 




0.9686 


0.7631 


Optimal encoding & recovery: 1 iteration 


T^al, C a i 


0.9719 


0.7805 


Optimal encoding & recovery: 100 iterations 




0.9997 


0.6261 


Optimal encoding, no recovery 


7^-0) Qa 


0.7445 


0.9091 


Optimal encoding & recovery: 1 iteration 


T^bl, Cbl 


0.7843 


0.9441 


Optimal encoding & recovery: 100 iterations 


T^&lOOj Cftioo 


0.7412 


0.9997 



Table 1 : Average fidelities 



For each of the optimizations, the process matrices Xc, Xr G C 8x8 , associated respectively 
with each C and 1Z were of reduced rank. For all the optimized C, each process matrix Xc was 
found to have a single dominant singular value, and hence, there is a single dominant 4x2 OSR 
element which characterizes C. For the optimized 1Z, each Xr was found to have two dominant 
singular values, and hence, there are two dominant 2x4 OSR elements which characterize 1Z. 
For example, the recovery /encoding pair (7£ Q i o, Caioo) has the OSR elements: 



-0.629 
0.455 + 0.378?' 
0.42 + 0.063?' 
0.13 + 0.233?' 



0.189 -0.332?' 
0.207 + 0.24?' 
-0.425 - 0.358?' 
0.626 + 0.226? 



Ri 



-0.707 0.532 - 0.342? 
0.134 + 0.087? 0.009 -0.166? 



0.194-0.175* 0.103 -0.138? 
-0.103 + 0.404i 0.833 - 0.276i 



-0.603 
0.313-0.103? 



-0.528 + 0.461? 
0.259 + 0.104? 



-0.262 - 0.131?' 0.172 -0.163? 
-0.374- 0.728? 0.174 - 0.333? 



It is not obvious that these correspond to any of the standard codes. However, by construction, 
C\Ci = I 2 and J2l=i — h- Referring to Figure |2j we can construct the encoding and 
recovery unitaries as, 

" Ri 



U c = [Ci C 2 



U 



R 



R-2 



where C 2 G C 2x2 is arbitrary as long as Uc is unitary, or equivalently, C\C 2 = and C 2 C 2 = 
I 2 . Observe that Ur is already a 4 x 4 unitary. 

Bar plots of the magnitude of the elements in the primal-dual pairs (Xc, Yc) and (Xr, Yr) 
corresponding to (7?. a i00) Caim) and (7?-6ioo, Cwoo) 316 shown in figures |4] and respectively. 
From many of such similar plots we have observed some common structure which may be used 
to reduce the computational burden. 
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Type 


K, C 




£ b 


Robust encoding, no recovery 


7^0, C a bl 


0.8840 


0.8840 


Robust encoding & recovery: 1 iteration 




0.9284 


0.9284 


Robust encoding & recovery: 100 iterations 


T^ablOO, C a b20 


0.9576 


0.9576 



Table 2: Average robust fidelities 



We also computed a robust encoding and recovery for the error set {£ a , £ } by iterating 
between (I2TT) and (l22b . The resultant average fidelities are in Table|2j 

Comparing Tables [l] and clearly shows that a robust design is possible although at a cost 
of performance. Also in this case after 100 iterations the robust fidelity did not increase. In ad- 
dition, the rank of the process matrices Xc and X R remained as before at 1 and 2, respectively, 
and the resulting OSR elements do not appear standard. 

6 Conclusions 

We have shown that the design of a quantum error correction system can be cast as a bi-convex 
iteration between encoding and recovery, each being a semidefinite program (SDP). We have 
also shown that the dual optimization, also an SDR is of lower complexity and thus requires 
less computational effort. The SDP formalism also allows for a robust design by enumerating 
constraints associated with different error models. We illustrated the approach with an example 
where the error system does not assume independent channels. 

Note added While this work was finalized for submission we became aware of the closely 
related [FSW06] and the subsequent commentary [RWA06]. 
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Figure 5: Magnitudes of primal-dual pairs (Xc, Yj) and (Xr, Yr) corresponding to 
(TZ bl00 , Cbioo) optimized for £ b . 
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A Dual problem 

We apply Lagrange Duality Theory [BV04, Ch.5]. Write the primal problem (l23l as a mini 
mization, 

minimize —TV XW 
Subject tO X > 0, Y^ij x ijCij = Im 

with optimization variable X = X^ G C nxn . The Lagrangian for (l27l) is, 



(27) 



L(X,Z,Y) = -TrXW = TrXZ-TrY(I m -^2 ij X ij BiB j ) 
= Ey Xvi-Wji - Z 3l + TV YCji) - TV Y 



(28) 



where Z = Z^ G c nXn and Y = Y^ G C mxm are Lagrange multipliers associated with the 
(Hermitian) inequality and equality constraints, respectively. The Lagrange dual function is 
then, 

g(Z,Y) = mi x L(X,Z,Y) 

-TV Y Zji — TV YCij — Wji (29) 
-oo otherwise 

For any Y and Z > 0, F) yields a lower bound on the optimal objective —TV X opt W. 
The largest lower bound from this dual function is then max { g(Z, Y) \ Z > }. Eliminating 
Z, this can be written equivalently as, 

minimize TV Y 

subject to KiY) - W > 0, K, tj (Y) = TV YCy { } 

with optimization variable Y = Y^ G C mxm . This is precisely the result in d24l . Because the 
problem is strictly convex, the dual optimal objective is equal to the primal optimal objective as 
stated in the first line of d25t . The complementary slackness condition gives the second line in 

d2J>. 
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